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Abstract 

In this paper, a linear univariate representation for the roots of a zero-dimensional polynomial 
equation system is presented, where the roots of the equation system are represented as linear 
combinations of roots of several univariate polynomial equations. The main advantage of this 
representation is that the precision of the roots can be easily controlled. In fact, based on the 
linear univariate representation, we can give the exact precisions needed for isolating the roots of 
the univariate equations in order to obtain the roots of the equation system to a given precision. 
As a consequence, a root isolation algorithm for a zero-dimensional polynomial equation system 
can be easily derived from its linear univariate representation. 

Key words: Zero-dimensional polynomial system, linear univariate representation, local 
generic position, root isolation 



1. Introduction 

Solving polynomial equation systems is a basic problem in the field of computational sci- 
ence and has important engineering applications. In most cases, we consider zero-dimensional 
polynomial systems. We will discuss how to solve this kind of systems in this paper. In par- 
ticular, we will consider how to isolate the complex roots for such a system. 

One of the basic methods to solve polynomial equation systems is based on the concept 
of separating elements, which can be traced back to Kronecker (114!) and has been studied 
extensively in the past twenty years (HI; [2; i; 0; i; [3; [ll|; [H; ES]^ H Ei)- The idea 



of the method is to introduce a new variable t = CiXi which is a linear combination of 
the variables to be solved such that t = CiXi takes different values when evaluated at 
different roots of the polynomial equation system 7^ = 0. In such a case, we say that t is a 
separating element for V ~ Q. U t ~ (^i^i is a separating clement for V ~ 0, the roots 
of P = have the following rational univariate representation (RUR): 

f{t) = 0,x,^R,it),i^l,...,n, 
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where / G Q[t] and Ri{t) are rational functions in t. As a consequence, solving multi-variate 
equation systems is reduced to solving a univariate equation f{t) ~ and to substituting 
the roots of f{t) =0 into rational functions Ri{t). Along this line, better complexity bounds 
and effective software packages for solving polynomial equations such as the Maple package 
RootFinding by Rouillier (|20l ) and the Magma package Kronecker by Giusti, Lecerf, and 
Salvy are given. 

The above approaches still have the following problem: for an isolation interval [a, h] of 
a real root a of f{u) = 0, to determine the isolation interval of x,; ~ Ri{a) under a given 
precision is not a trivial task. In this paper, we propose a new representation for the roots 
of a polynomial system which will remedy this drawback. 



Fig. 1. The distribution of the roots of Ti{x) = 0(i = 1, 2, 3). The red diamonds (blue 
crosses, black circles) are roots of Ti{x) = (T2{x) — 0, T3{x) = 0) and red (blue) 
boxes are neighborhoods for the red diamonds (blue crosses). 

In the ISSAC paper (0), based on ideas similar to separating elements, a local generic 
position method is introduced to solve bivariate polynomial systems and experimental re- 
sults show that the method is quite efficient for solving equation systems with multiple 
roots. In this paper, we extend the method to solve general zero-dimensional polynomial 
systems. A local generic position for a polynomial equation system "P = is also a linear 
combination of the variables to be solved: t = CiXi which satisfies two conditions. First, 
tk = X^iLi '^i^i is a separating element of Vk = {V) D Q[xi, . . . , Xk] for k = 2, . . . , n, and 
the roots of Vk = have a one-to-one correspondence with the roots of a univariate equa- 
tion Tk{tk) = 0. Second, for a root ^ = (^i, . . . of Vk = represented by a root i] of 
Tk{tk) — 0, all the roots rjj of Tk+i{tk+i) = corresponding to the roots of Vk+i ~ 0, say 

~ (^,^fc+ij), "lifted" from ^ are projected into a fixed square neighborhood of r\. This 
"local" property is illustrated in Figure 1. We prove that if t.^ = X]"=i "^i^j ^ local generic 
position for 'P, then the roots of P = can be be represented as special linear combinations 
of the roots of univariate equations Tk{tk) = 0, fc = 1, . . . , n: 

a2-ai a„-a„_i 
|(ai, , . . . , ) I Tk(ak) = 0}, 

Sl Sl---S„_l 



where Sj are certain positive rational numbers and the ctj+i matching aj are in certain 
square neighborhood of aj to be defined in Section 2. Such a representation is called a 
linear univariate representation (LUR for short) of the polynomial system. 

The main advantage of the LUR representation is that the precision of the roots can be 
easily controlled. For RUR, computing solutions with a given precision is not a trivial task 
as we mentioned before. It is not easy to know with which precision to isolate the roots 
of f{t) = is enough in order for the roots of the system Xi = Ri{t) to satisfy a given 
precision. For LUR, precision control becomes very easy. We can give an explicit formula 
for the precision of the roots of Ti (x) = in order to obtain the roots of the system with 
a given precision. So we can obtain the solutions of the system by refining the roots of 
Ti{x) = at most once. Another advantage of LUR is that when we isolate the roots of 
Ti+i{x) = 0, we need only to consider a fixed neighborhood of each root of Ti{x) = 0. 

We propose an algorithm to compute an LUR for a zero-dimensional polynomial system. 
The key ingredients of the algorithm are to estimate the root bounds oi V ~ and to 
estimate the separating bounds for the roots of Vk+i = lifted from a root of Vk = 0. We 
adopt a computational approach to estimate such bounds in order to obtain tight bound 
values. For the root bounds of T-" = 0, we use Grobner basis computation to obtain the 
generating polynomial of the principal ideal (V) riQ[xi\ and use this polynomial to estimate 
the root bound for the Xi coordinates of the roots of = 0. The separating bounds for 
Vk ^ are obtained from the isolating boxes for the roots of the Tk{x) = 0. These bounds 
in turn will be used to compute the isolating boxes for the roots of Vk+i = 0. Hence, the 
algorithm to compute an LUR also gives a set of isolating boxes for the roots of = 0. 

In Section 2, we give the definition of LUR and the main result of the paper. In Section 
3, we present an algorithm to compute an LUR of a zero-dimensional polynomial system as 
well as a set of isolation boxes of the roots of the equation system. In Section 4, we provide 
some illustrative examples. We conclude the paper in Section 5. 

2. Linear univariate representation 

In this section, we will define LUR and prove its main properties. Let 

V = {fl{xi, . . . ,Xn), ■ . ■ , fs{xi, . . . , Xn)} 

be a zero-dimensional polynomial system in Q[.ti, . . . ,a;„], where Q is the field of rational 
numbers. Let 

I, = {V)nQ[xi,...,x,],i^ 
where {V) is the ideal generated by V. We use VciV) to denote its complex roots in C". 

Since we will use rectangles to isolate complex numbers, we adopt the following norm 
for a complex number c = x + yi: 

\c\ = max{|a;|, \y\}. (1) 

The "distancJ*}' between two complex numbers ci and C2 is defined to be |ci — C2|. It is easy 
to check that this is indeed a distance satisfying the inequality |ci — C2 1 < |ci — C3 1 -I- |c3 — C2 
for any complex number C3. Let cq be a complex number and r a positive rational number. 



* The results in this section are also valid if we use the usual distance for complex numbers. 



Then the set of points having distance less than r with cq, that is {ci S C | |ci — co| < r}, 
is an open square with cq as the center. 
By an LUR, we mean a set like 

{Ti{x), . . . ,Tn{x),Si,di,i = 1}, (2) 

where Ti{x) £ Q[a;] are univariate polynomials, s,; and di are positive rational numbers. The 
roots of (2) are defined to be 

, n and 

(3) 

where sq = 1- Geometrically, we match a root of r,;(2:) ~ with those roots oiTi^i{x) = 
inside a squared neighborhood centered at ai. See Figure 1 for an illustration. An LUR 
for 7-" is a set of form (2) whose roots arc exactly the roots oiV = Q. 

It is clear that an LUR represents the roots of V as linear combinations of the roots of 
some univariate polynomial equations. The LUR representation has the following advantage: 
wc can easily derive the precision of the roots of = from that of the univariate equations 
as shown by the following lemma. 

Lemma 1. Let (2) he an LUR for a polynomial system V ~ 0. Lf ai is a root of Ti{x) = 
0(1 < i < n) andlxi is an approximation of ai with precision e^, then the approximate root 
(ai, . . . , ofV = Q has precision max{ei, . . . , 

Proof. Since Xi = "'"i ^'^"^ approximate root of ai has precision e^, the approx- 
imate root Ti = has precision no larger than '^^^'^^-'^ ^ ■ 

For a zero-dimensional polynomial system V, let di,ri (i — and s; (i — 

1, . . . , 71 — 1) be positive rational numbers satisfying 

A = min{^|a ~ Vr; G Vc{I^^,), (r;, a), (ry, /3) £ Vc{I^), a ^ /3}, 

< min{A, 4^-^}, 

ri > 2max{|ai|,V(ai, . . . ,ai) g Vc{Tt)}, 
di 

St < 

ri+i 

where sq = 1, do = +oo. Geometrically, Di is half of the root separation bound for roots of 
Xi considered as points on a "fiber" over each root of is twice the root bound for 

the i-th coordinates of the roots of li, and Si, the inverse of the slope of certain line, is a 
key parameter to be used in our method. If V?7 G Vc{Xi-i), #{a|(r/,a) € Vc{Ii)} = 1, we 
can choose any positive number as di. 



a2 — ai 

Sl 



Sl • • ■ Sn-1 
•Sj_idj,i = 1 



)|T,(aO = 0,i = l,, 
,n-l} 



(4) 

(5) 
(6) 
(7) 



For Si satisfying (7), consider the ideal 

ii = (li U {x - xi - S1X2 si • • • Si-ixJ), (8) 

where a; is a new variable. It is clear that 1; is a zero-dimensional ideal in Q[xi, . . . ,Xi,x]. 
And the elimination ideal (li) D Q[x] is principal. Let Ti{x) be the generator of this ideal: 

(lOnQM - (r,(x)). (9) 

The following is the main result of this paper. 

Theorem 2. If di,Si satisfy conditions (5), (7) and Ti is defined in (9), then the corre- 
sponding set (2) is an LUR for V . 

We will prove two lemmas which will lead to a proof for the theorem. For a root r]i of 
Ti{x) = 0, let 

S,,, = {7? e C I I77 - 77i| < pi},i = l,...,n (10) 

where pi = si • • • Si-idi, (sg = !)■ Note that S,,. is an open square whose center is rji and 
whose edge has length 2/9,;. With this notation, the roots of (2) can be written as 

{[ai, , . . . , ) I J,(q;,) = 0,1 = 1, . . . , n and 

Si Si • • • Sn~l 

Qfi+i G Sa;, i = 1, . . . , n - 1} (11) 
In Figure 1, are interior parts of the squares. We have 

Lemma 3. Under assumptions of Theorem 2, we have S'rii+i C Sr^^i^l,- ■ ■ ,n-l^ where 
(Ci, ■ • ■ , ii+i) e Vcili+i) and 

'7j==Ci + si6H h si • • • s,_i (12) 

Vi+i = 6 + S16 H h si • ■ • s,_i + si • ■ • s,6+i = ?/j + si • ■ • s,^i+i. (13) 

Proof. From the definition of in (8), rji is a root of Ti{x) ~ 0, is a root of Ti^i{x) = 0, 
and each root of Ti^i{x) = has the form (13). 

We first prove that rji+i G Using (6) and (7), we have 

-Vi\= si- \ < ^si • ■ -Siri+i < ifii • --Si-idi = ^pt. (14) 

As a consequence, iji^i is in S,,. . 

We now prove that S,,j_|_j C Sr^ - By (5), we have p^+i = si • • • s^iii+i < ^si ■ • • Si^idi — 
\pi. Therefore, for any 77 e by (14), we have I77 - ■qi\ < \r] - + - < 

Pi+i + ^Pi < Pi- Hence rj G and the lemma is proved. | 

For rational numbers Oj , we call fi = E*-^ 2 CLj X 7" Qb separating element of 2i, if Va, /3 G 
Vcili), 13 implies fi{a) ^ fi{j3) (see paper (|20)). 
Theorem 2 follows from (d) of the following lemma. 



Lemma 4. Under assumptions of Theorem 2, for i ~ 1, . . . ,n, we have 
(a) X ~ xi + si X2 + ■ ■ ■ + si ■ ■ ■ Si-iXi is a separating element ofli. 

(h) Each root rji of Ti{x) = is in an for a root rji-i of Ti^i{x) — 0. Furthermore, 

if Vi-i = ^1 + si ^2 + ■ ■ ■ + Si ■ ■ ■ Si-2^i-i. then all roots of Ti{x) — in ^rn-i Oife of the 
following form 

Vi Vt-i + si • • ■ Si-i (15) 

where (^i, . . . , ^j) G Vc{Ii). 

(c) Srii are disjoint for all roots rji of Ti{x) = 0. 

(d) (Ti (x), . . . ,Ti{x),Sj,dj,j ~ 1, . . . ,i ~ 1) is an LUR for li . 

Proof. We will prove the lemma by induction on k = i. For k = 1, since (Ii) = (Ti(a;)), 
statements (a) and (d) are obviously true. We do not need prove (b). From (5), we have 
di < min{i|Q! — /3|,Va,/3 € Vc(Ii) = Vc(Ti),a /?}. As a consequence, S^^ are disjoint 
for all roots rji of Ti(x) = 0. Statement (c) is proved. 

Suppose that the result is correct for fc = 1, . . . , i. We will prove the result for fc = i + 1. 

We first prove statement (a). Let ^ = (^i, . . . , Ci+i) a-nd /3 — . . . , A+i) be two distinct 
elements in Vc{T^i+i)- We consider two cases. If (^i,...,^;) is different from (/3i, . . . , /3i), 
then by the induction hypothesis rji = + si^2 + • • ■ + si ■ ■ ■ is also different from 

(^i = Pi+siP2 + - ■ • ■ ■ Si_i/3i. By (c) of the induction hypothesis, S^^ and Sg. are disjoint. 
By Lemma 3, rj^+i = t^j + si • • ■ Si^i+i £ S^, and 9i+i = 9i + si - ■■ Si(3^+i G Sg^. Then, in 
this case we have rji^i ^ 6*^+1. In the second case, we have = (/3i, . . . , /?;). 

Then, r]i = 9i and ^i+i ^ Pi+i- It is clear that 77^+1 = ?7i + si • • • Si^i^i is different from 
di+i = 9i + si ■ ■ ■ Sif3i+i. Thus, (a) is proved. 

We now prove statement (b). Use notations in (12) and (13). By Lemma 3, we have 
?7i+i G Then, each root of Ti^i{x) = is in an for a root 77, of Ti{x) = 0. Let 
(/?i, . . . , Pi+i) G Vc(Ii+i) such that 9i^i = /3i+si/32 + - ■ - + 51 ■ • • Si/Si^i is another clement in 

. We claim that (/3i , . . . , /3i) must be the same as (^1 , ... , ^i). Otherwise, by the induction 
hypothesis (a), 9i ~ (3i + si/32 + • ■ • + si ■ ■ • si^iPi is different from rji. By the induction 
hypothesis (c), S^. and Sg. are disjoint which is impossible since by Lemma 3, 9i^i G 
and 9i+i G Thus, (/3i, . . . ,^i) = (^1, . . . ,^j) and hence 9i+i = 77^ + si • ■ • s^Pi+i. This 
proves equation (15) and hence statement (b). 

We now prove statement (c). Use notations in (12) and (13). By Lemma 3, ^rn+i C S,,. . 
As a consequence, we need only to prove that the squares Srji+i contained in the same Sri^ 
arc disjoint. Let ?7i+i, be two roots of Ti+i{x) = in . By statement (b) just proved, 
we have 

Vi+i =m + Si - ■ ■ s^^i+i, 6*4+1 ='ni + Si - ■ ■ Sij3.i+i 

where 77,; is defined in (12) and (^1, . . . ,^1,^^+1), (^1, . . . are roots of Ii+i. Then, 

by (5), 

hi+i - 9i+i \ = si • • • Sj|Ci+i - Pi+i \ >2si - ■ ■ Si dj+i = 2pi+i. 
So, §r)i+i and ^Oi^i are disjoint. Statement (c) is proved. 

Finally, we prove statement (d). Let ^ = (^i, . . . ,Ci+i) € Vc(2Ii+i) and 77^ = ^1 + si^2 + 
• • • + si • • • Sj-i = 1, . . . , z + 1. By the induction hypothesis, we have (^1, . . . S,i) — 
('?!' ^^^IT^' • ■ • ' s'i~X~-i ) "^^^^^ - V]\ < si • • • Sj-idj,j = 1, . . . , i. Note that the in- 

equality is equivalent to that rjj^i G S^^ • By (15), we can recover the ^i+i with the following 
equation 

^ 'n^+l - V, 



From Lemma 3, we have rji^i G E,-/- or cquivalently \r]i+i — ?7i| < si • • • Si-idi. Then the root 



(6, • • •e.+i) = (^71, ■ • • , '-^j^) is a root of the LUR {T^{x), . . . , r,+i(x), s„ d„ j = 



!,...,«). We thus proved that the roots of li+i are the same as the roots of the LUR and 



We have the following corollaries. 

Corollary 5. If (2) is an LUR for a polynomial system V , then the roots of Xi — are in 
a one to one correspondence with the roots of Ti{x) = for i ~ 1, . . . ,n. 

Proof Let ^ = (^1, . . . , ^i) e Vc{Ii). Then iji = + si^2 H h si • • • Si-i is a root of 

Ti{x) = 0. By (a) of Lemma 4, this mapping is injective. This mapping is clearly surjective. 
I 

Corollary 6. The real roots ofV~0 are in a one to one correspondence with the real roots 
of Tn{x) = 0. More precisely, if an is a real root of Tn{x) = 0, then in the corresponding 



root {ai, "^^"^ , ■ ■ ■ , ^ ) ofV = 0, ai is a real root ofTi{x) = 0, i = 1, . . . , n — 1. 



Proof. For each root r] of Ti_i (cc) = 0, let be the open square neighborhood of 77 defined in 
(10). We claim that a real root of Ti{x) = cannot be in S,, for a complex root rj of Ti_i(x) = 
0. Since Ti-i{x) has rational numbers as coefficients, the complex roots of Ti^i{x) ~ 
appear as pairs which are symmetric with the real axis and the open square neighborhoods 
for a pair of complex roots are disjoint. Then the open square neighborhood of any complex 
root has no intersection with the real axis. This proves the claim. As a consequence, if a„ 
is a real root of T„(a;) — 0, then a„ is in the open square neighborhood of a real root Q!„_i 
of Tn-i{x) = 0. Repeating the process, we obtain a real root (ai, , . . . , ) for 

V ~ where all ai are real numbers. The other side is obvious: a real root of V ~ will 
correspond to a real root of Tn{x) =0. | 

From the lemma, we can consider the real roots of an LUR if we only interest in the real 
roots of T' = 0. 

3. Algorithm for computing an LUR and roots isolation 

In this section, we will present an algorithm to compute an LUR for a zero-dimensional 
polynomial system. The algorithm will isolate the roots of the system in C" at the same 
time. 

3.1. Complex isolation intervals and isolation boxes 

We introduce some basic concepts of interval computation. For more details, we refer to 



Let DQ denote the set of intervals of the form [a,b], where a < b <E Q. The length of 
an interval / = [a, b] € UQ is defined to be |/| = b — a. Assuming ai < 02, we define the 
distance between two intervals as 



hence statement (d). 



I 




A pair of intervals (/, J) is called a complex interval, which represents a rectangle in 
the complex plane. A complex number (a, p) = a + /3i (i^ = —1) is said to be in a complex 
interval (/, J) if a G / and 13 ^ J. The length of a complex interval (/, J) is defined to be 
|(/, J)| = max{|/|, I J|}. We define the distance between two complex intervals as 

Dis(([ai,6i], [pi,qi]), ([a2,&2], [P2,q2])) = max{Dis([ai, 61], [02, b2]), Dis([pi, gi], [^2,92]}. (16) 



A set S of disjoint complex intervals is called isolation intervals of T{x) = if each 
interval in iS contains only one root of T{x) =^ and each root of T{x) = is contained in 
one interval in S. Methods to isolate the complex roots of a univariate polynomial equation 
are given in (0; [ll [H [ll) . 



Let DC denote the set of complex intervals. An element {lf,ll) x ■ • • x (I^Ji) in DC" 
is called a complex box. A set S of isolation boxes for a zero dimensional polynomial 
system V in Q[a;i, . . . , Xn] is a set of disjoint complex boxes in DC" such that each box in S 
contains only one root of = and each root of "P = is in one of the boxes. Furthermore, 
if each box B = (If, /J) x • • • x (I^, l\) in S satisfies max{|/f |, < e, then S is called an 

i 

e-isolation boxes of 'P = 0. The aim of this paper is to compute a set of e-isolation boxes 
for a zero-dimensional polynomial system V. 

3.2. Grobner basis and computation of ri and Ti{x) 

In this subsection, we will show how to use Grobner basis to compute defined in (6) 
and Ti{x) defined in (8) supposing the parameters Si are given. 

Let 7-* C Q[xi, . . . , Xn] be a zero-dimensional polynomial system. Then A = Q[xi, . . . , Xn]/ 
(V) is a finite dimensional linear space over Q. Let Q he a Grobner basis of V with any 
ordering. Then the set of remainder monomials 

B = {x*^ • • • x*n \x*i ■ • ■ x*n is not divisible by the leading term of any element of Q} 

forms a basis of ,4 as a linear space over Q, where ti are non-negative integers. 
Let / G . . . , Xn]. Then / gives a multiplication map 

Mf.A — >A 

defined by Mf{p) = fp for p Q A. It is clear that Mf is a linear map. We can construct the 
matrix representation for Mf from B and Q. The following theorem is a basic property for 

Mf dil). 

Theorem 7 (Stickelberger's Theorem). Assume that V C Q[a;i, . . . , x„] has a finite positive 
number of solutions over C. The eigenvalues of Mf are the values of f at the roots ofV = 
over C with respect to multiplicities of the roots ofV = 0. 

Let Si be rational numbers satisfying (7) and 

^ V U {x ~ xi ~ S1X2 - . ■ • - si • ■ • Si-iXi}. 

We can compute gi{xi) and Ti{x) such that 

{9^{xi)) = Q[x,] n {V) and {T,{x)) ^ Q[x] n [T,). (17) 

In fact, we can construct the matrixes for M^,. and M^ based on B and Q, and gjixi) and 
Ti(x) are the minimal polynomials for M^^ and Mj,, respectively (See reference (|6|)). Note 
that we can also use the method introduced in reference (0) to compute gi{xi),Ti{x). 



From Theorem 7 and (a) of Lemma 4, the i-th coordinates of ah the roots of 7^ = are 
roots of gi{xi) = 0, and all the possible values oi x = Sj=i si ■ • • Sj-iXj on the roots of 
7^ = are roots of Ti{x) = 0. 

Now we show how to estimate defined in (6). At first, compute {gi{xi)) — {V) nQ[xi]. 
Then we have the following result. 

Lemma 8. Use the notations introduced before. Then 

= 2max{RB(5,(x,))} (18) 

satisfies the condition (6), where RB((7) is the root bound of a univariate polynomial equation 
9 = 0. 

Proof. The lemma is obvious since for any root {£,1 , . ■ . , € Vc{Ti) , £i is a, root of gi{xi) = 



3.3. Theoretical preparations for the algorithm 

In this subsection, we will outline an algorithm to compute an LUR for V and to isolate 
the roots of T-* = under a given precision e. The algorithm is based on an interval version 
of Theorem 2. 

We define the isolation boxes for an LUR defined in (2) as: 

{Si X X ■ ■ ■ X ^" ~ I B, e S,, Dis(B,+i, B,) < p,/2, 1 < i < n - 1} (19) 

Si Si • ■ • S„_i 

where Bi is a set of isolation boxes for the complex roots of Tiix) — and = si ■ • • si^idi. 
In Theorem 16 to be proved below, we will give criteria under which the isolation boxes for 
V are the isolation boxes of an LUR. 

Let V C Q[a:i, . . . , a;„] be a zero-dimensional polynomial system. We will compute an 
LUR for V and a set of e-isolation boxes for the roots of P = inductively. 

At first, consider i — 1. We compute Ti{x) as defined in equation (17). Let Bi be a set 
of isolation intervals for the complex roots of Ti{x) = 0. Then, we can set di to be the 
minimal distance between any two intervals in Bi. 

For i from 1 to n — 1, assuming that we have computed 

• An LUR (Ti(a;), . . . , Ti(a;), Sj, dj, j = 1, . . . , i — 1) for 2^. 

• A set of e-isolation boxes for Xi . 

• The parameter di. 

We will show how to compute ri+i, s^, TiJ^iix)^ dj+i, and a set of e-isolation boxes of 
the roots of Ii+i = 0. The procedure consists of three steps. 

Step 1. We will compute r^+i, Si as introduced in (6) and (7). With Si, we can compute 
Ti^iix) as defined in (17). 

Here ri+i can be computed with the method in Lemma 8. Note that di is known from 
the induction hypotheses. Then we can choose a rational number Sj such that condition 
(7) is valid. Finally, Ti+i[x) can be computed with the methods mentioned below equation 
(17). 



Step 2. We are going to compute the isolation intervals of the roots of li+i = 0. Let 
^ = (^1, . . . , ^i) be a root of 1; = 0. We are going to find the roots of li+i = "lifted" from 
^, that is, roots of the form 

= (Ci, ■ • • ,6,6+1 j), J = 1,. . .,m. (20) 

To do that, we need only to find a set of isolation intervals for 6+i,j with lengths no larger 
than e, since we already have an e-box for ^. 
Let 

?7i = + si6 H h si • • • Si-i6- 

Then, jy,; is a root of Ti{x) ~ 0. By (b) of Lemma 4 the roots of Ti_|_i(a;) = correspond 
to C,j are 

0j=r]i + si--- Siii+ijJ 1, . . . ,m. (21) 
We have 

Lemma 9. Let Ii = ([a, 6], [c, d]) 6e an isolation interval for the root rji of Ti{x) = such 
that \Ii\ < jpi where pi = si ■ ■ ■ Si-idi. Then all 9j in (21) are in the following complex 
interval 

I/, =((a-p,/2,6 + p,/2,(c-p,/2,d + p,;/2)). (22) 
Furthermore, the intervals are disjoint for all roots i] of Ti{x) = 0. 

Proof. In Figure 2, let square ABCD be — {6 e C| |0 — ?7i| < pi} and square 
AiBiCiDi = {6* e C I 16* - ?7i| < pj/2}. By equation (14), we know \0j - < \pi. So, 
Oj is inside AiBiCiDi. Let rectangle A2B2C2D2 be the complex interval and rectangle 
A3B3C3D3 the complex interval I/, which is obtained by adding pi/2 in each direction of 
the rectangle A^B^C^D^. Then, I/, contains AiBiCiDi and hence 9j is inside I/^. 

For any 6* S I/. , we have |6' — ry^l < \9 — P\ where P is one of the points ^2, -^2, C2, -D2 to 
make 16* -P| maximal. It is clear that \9-P\< pi/2+\Ii\ = fp^. So, V6I G I/, , |6I - r?i| < 
Since is the set of complex numbers satisfying \9 — r]i\ < pi, we have I/^ C S^^. By (c) 
of Lemma 4, S^^ arc disjoint for all roots of Ti{x) = 0. Then Ij. arc disjoint for all roots of 
Ti{x) ^ too. I 

The following lemma shows what is the precision needed to isolate the roots of Ti+i(a;) = 
in order for the isolation boxes to be contained in some I/. . It can be seen as an effective 
version of the fact r\ij^x S S^. proved in Lemma 3. 

Lemma 10. Use the notations introduced in Lemma 9. Let {Bj, j ~ 1, • • . , rii} be a set of 
^Pi-isolation boxes for the roots 9j,j = 1, . . . , m of Ti^i(x) = 0. Then, for all j 



B-j C I/, and Dis{Bj,Ii) < pi/2. 



(23) 
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Fig. 2. The isolation intervals /i, S^. , I/, for a root rji of Ti(x) — 0. 
r]i is represented by o. 

Proof. From the proof of Lemma 9, the distance from 77^ to line BC in Figure 2 is pi and 
the distance from r\i to line Bt,Cz is less than So, the distance between line BC and 
BaCa is at least -jpi. This fact is also valid for the pairs of lines AD / A^^D^. AB/A3B3, and 
CD/C3D3. Since the isolation boxes Bj are of size smaller than pi/4 and their centers are 
inside A3B3C3D3, the boxes Bj must be inside ABCD. Note that li is rectangle A2B2C2D2- 
Since 9j is inside both and rectangle A^iB^CzD^ and the distance from r/^ to each edge 
of A^B^C^Dj, is pi/2, the distance between Bj and li must be smaller than pi/2. | 

Isolate the roots of Ti+i(a;) = with precision jpi. By Lemma 10, all the roots are in 
one of the intervals 1/ where / is an isolation interval for a root 77 of Ti(yX) = 0. 

Let Kj = {[pj^qj],[gj,hj]){l < j < m) be the isolation intervals for the roots Oj of 
Ti+i{x) = inside I/^. Then from (21), the isolation intervals of <^i+ij(l < j < m) are 

j_ ^ . ^ - h ^ {[pj -b,qj -a],[gj -d,hj - c]) 
We have 

Lemma 11. With the notations introduced above, if the following conditions 

{qj ~ Pj) + {b - a) < si ■ ■ ■ Sie, {hj - gj) + {d - c) < Si ■ ■ ■ s,e (25) 

Srj. = ^pin Dis(([pfc, [gt, hk]), {[Pj,qj], [9j,hj])) > ma.x{b -a,d- c}. (26) 
are valid, then Ji+ij defined in (24) are e-isolation intervals o/^j+ij- in equation (20). 



Proof. It is clear that condition (25) is used to ensure the precision: j | < e. 

We consider (26) below. Assume that Ji+i,fc(l < k j < m) are any two intervals 

defined in (24) for the (i + l)-th coordinates of the roots (^1, . . . , £_i,£,i^ij), (^1, . . . , Cii Ci+i,/c) 
of li+i = 0, respectively. Since we have derived the e-isolation boxes for the roots of 
li = 0, we need only to ensure that the intervals of the {i + l)-th coordinates of the roots 
of = lifted from a fixed root of 1; =0 are isolation intervals. That is, to show 
Dis(Ji+ij, Ji+i,fe) > 0. 

Assume that Kj = {[pj,qj],[gj,hj]) and Kk = {[Pk,<}k],[gk,hk]) are the isolation in- 
tervals of the roots rjj, rjk of Ti+i{x) = 0. Here rjj, rjk correspond to (^1, . . . , ^i+i^). 



(^1, . . . , ^i, respectively. So Kj, Kk correspond to Ji+ij, >/i+i,fe, respectively. Assume 

that pj < pk , gj < gk ■ Then we have 

max{Dis([pj - b, qj - a], [pk - b, qu - a]), Dis([c;j - d, hj - c], [gk - d, hk - c])} 

UlS^Ji+lj-, Ji+l.fcJ = , 

Sl - ■ ■ Si 

and 



Ci =Dis{[pj-b,qj-a],[pk-b,qk-a]) = 
£2 = Dis([g_,--d, hj-c], [gk~d, hk-c]) = ■ 



[Pk - gj) - (fe - a), if (Pfc - qj) - (& - a) > 0, 
0, otherwise, 



I {9k - hj) -{d- c), if {gk - hj) - (d - c) > 0, 
I0, otherwise. 
Kj and are disjoint since they are isolation intervals of T'i+i(a;) = 0. So 

I)is{Kj,Kk) = max{pfc - 17^,3^ - /ij} > 0. 

It is clear that Dis(Ji+i.j, Ji+i.k) > if £1 > or £2 > 0. Then we conclude if inequality 
(26) is true, then Dis( J^+ij, Ji+i.k) > 0. This proves the lemma. | 



Geometrically, 5^. is the separation bound for the roots of Ti+i(x) = corresponds to 
those roots of li+i lifted from the root of 1^ = corresponding to the root rji of Ti{x) = 0. 

Remark 12. Note that in (26), we obtain /.j = ([a, b], [c, d]) first and Kj — {[pj, qj], [gj, hj]) 
later. We will refine the isolation interval li of Ti{x) = such that inequality (26) is true. 
After the refinement, all other conditions are still valid. We need to do this kind of refinement 
only once. 

As a consequence of the above lemma, we have 

Corollary 13. Let B be an e-isolation box for the root ^ o/I^ = and Ji+i.j defined in 
(24)- If (25), (26) are valid, then B x Ji+ij,j — 1, . . . , m are e-isolation boxes for the roots 
Pj o/Xi+i = 0, which are lifted from ^. 

Step 3. We will show how to compute di+i from the isolation intervals of Ti+i{x) = 0. 

Lemma 14. Let ^ 

min{-— ^i^, -^}, (27) 
2si • • • Si ISi 

where Si+i is the minimal distance between any two isolation intervals ofTi+i{x) — 0. Then 
di+i satisfies conditions (5). 



Proof. Let 6j and Ok be two different roots of Ti+i{x) = defined in (21). Then we have 

^ ^ _ Oj - Ok 

~ <ii+l.k — ■ 

Si...Sj 

Therefore, -Di+i — imnri£Vc{Ti{x}){ 2s^^- -s ■ i ^^'^ parameter defined in (4), where Srj is 
determined as in (26). It is clear that Di+i is not larger than Si+i which is the minimal 
distance between any two isolation intervals of Tj+i(x) = 0. Then, the first condition in 
(5) is satisfied. In order for the second condition in (5) to be satisfied, we also require 
di+i < So the lemma is proved. | 



We can summarize the result as the following theorem which is an interval version of 
Theorem 2. 

Theorem 15. Let (2) be an LUR such that di, Vi, and Si satisfy (27), (6), and (7) respec- 
tively, Bi the ei-isolation boxes for the roots ofTi(x) — 0, and Si = min{Dis(_Bi, B2) \ Bi, B2 G 
B,,B,^B2}. If 

ei < e,ei + e.,+i < si ■ ■ ■ s^e, < ^, e,+i < ^, < S'i+i, (28) 
where Pi = si ■ ■ ■ Si-i di, then (19) is a set of e-isolation boxes for V ■ 

Proof. We first explain the functions of the inequalities in (28). The first two inequalities 
in (28) are introduced in (25) to ensure the e precision for the isolation boxes. The third 
inequality in (28) is introduced in Lemma 9 to ensure 9j £ I/, and li- are disjoint. The 
fourth inequality is introduced in Lemma 10 to ensure the isolation intervals of the roots of 
Ti^i{x) = are inside their corresponding interval I/^. The last inequality is introduced in 
(26) to ensure the recovered isolation boxes of V are disjoint. 

Now the theorem is a consequence of Corollary 13. Here, we give the explicit expression 
for the isolation boxes. The expression for interval Ji+ij in (24) is directly given. The 
matching condition Dis{Bi+i, Bi) < pi/2 is from condition (23). | 

We have the following effective version of Theorems 2 and 15 by giving an explicit formula 
for e,;. 

Theorem 16. Using the same notations as Theorem 15. Let e be the given precision to 
isolate the roots of V . Let 

■ f sie di 
ei =mm{e, — , ^,'^2), 

. rSi---Si-ie si---Sie Si ■ ■ ■ s^-idi Si • • • Si_2di_i , , , 

=mm{ , , , ,S^+i}, (29) 

where i = 2, ...,n, sq = 1,s„ = 1, Sn+i = +00. If we isolate the roots of Ti{x) ~ with 
precision ei, then (19) is a set of e-isolation boxes for V = 0. 

Proof. By (29), we have e,; < ^J-^^ and e^+i < ^J-^^. Then the second inequality in (28), 
Ei + Ei+i < si • • • SiC, is valid. All other inequalities in (28) are clearly satisfied and the 
theorem is proved. | 

We can also compute the multiplicities of the roots with the LUR algorithm. 

Corollary 17. If we compute the last univariate polynomial Tn{x) in the LUR as the char- 
acteristic polynomial of Mx, then the multiplicities of the roots ofP^O are the multiplicities 
of the corresponding roots of Tn{x) = 0. 

Proof. By (a) of Lemma 4, 2; = xi + si X2 + • ■ • + si • • • s„_ia;„ is a separating element. By 
Theorem 7, the characteristic polynomial of keeps the multiplicities of the roots of the 
system. The corollary is proved. | 



3.4- Algorithm 



Now, we can give the full algorithm based on LUR. 

Algorithm 1. The input is a zero dimensional polynomial system V = {Pi,...,Pt} in 
Q[xi, . . . , Xn] and a positive rational number e. The output is an LUR for V and a set of 
e-isolation boxes for the roots of "P = 0. 

51 Compute a Grobner basis Q oi V with the graded reverse lexicographic order and a 
monomial basis B for linear space A = Q[xi, . . . ,Xn]/ {V) over Q. 

52 Compute Ti{x) as defined in (17) with the method given in Section 3.2; compute a set 
of e-isolation boxes Bi for the complex roots of Ti{x) = 0; set di = min{Dis(i3i, B2) \ Bi, 

53 For i = 1, . . . , n — 1, do steps S4-S9; output the set of boxes (19). 

54 Compute Vi^i with the method in Lemma 8. Select a rational number Si such that 
condition (7) is valid. 

55 Compute Ti^i(x) as defined in (17) with the method given in Section 3.2. 

56 Set Pi = si • • ■ Si^idi and compute a set of ipi-isolation boxes Bi-^i for the complex 
roots of Tij^i{x) ~ 

57 Set S.,+i = min{Dis(Si, B2) | -Bi, B2 G B^+i,Bi ^ B2}. 

58 Compute d^+i with formula (27). 

59 Compute with formula (29); refine the isolation boxes Bi of Ti{x) = with precision 
ei; still denote the isolation boxes as Bi. 

Remark 18. From Lemma 9, the roots of Ti+i{x) — arc in the rectangle Ii- . So, we need 
only to isolate the roots of Ti(x) ~ inside these rectangles. This property is very useful in 
practice, see Figure 1 for an illustration. 

4. Examples 

In this section, we will give some examples to illustrate our method. 

We first use the following example to show how to isolate the roots of a system with our 
method. Note that with an LUR, we can also use floating point numbers to compute the 
roots of 7-" = if the floating point number can provide the required precision as shown in 
the following example. 

Example 19. Consider the system P := [x^ + i/^ + — 3, + 2y^ — 3z + l,x + y — z]. 
The coordinate order is {x,y, z). 

The Grobner basis Q with the graded reverse lexicographic order z > y > x of P is: 

[-x-y + z,x^ +2yx + 3x~4: + 3y,~3x + x^ + l-3y + 2y^,6x^ -30 + 9x^ + 25y + 5x]. 

The leading monomials of the basis are {z,xy,y^,x^}. So the monomial basis of the quotient 
algebra A = Q[xi, ...,Xn]/{P) is B = [l,x,y,x'^]. 



Let t = x, we can compute: 
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The minimal polynomial of Mt is 

Ti{t) = 5 - 60t + 6t^ + 18t^ + 6t^. 

Compute its complex roots with the function "Analytic" in Maple package [RootFinding] , 
we obtain 

Ri = [-2.22081423399575- 1.535197796461521,-2.22081423399575 
+ 1.535197796461521,0.0842270424726020,1.35740142551890]. 

Computing the roots distance with formula (16), we obtain di < 0.6365871918. We can set 

In a similar way, we compute My and its minimal polynomial g2iy) = —29 — 66'!/ + 
60 + 12 y^. The root bound of .92(2/) is 3. So we have r2 = 6. Since ^ = we set 

Let t ^ X + siy. We can compute a matrix Mt and its minimal polynomial 
T2{t) = 863337 - 6119640t+ 360000 + 1920000 + 640000 
Computing its complex roots, we have 



i?2 = [-2.24194942371773 - 1.413423955527621, -2.24194942371773 
+ 1.413423955527621, 0.143249906267126, 1.34064894116850]. 

Computing the minimal distance between any two roots, we have 5*2 = 0.5986995174. From 
equation (27), we can obtain 

, . r S2 di 
d2 = mm{ — , — } = 5. 
2 si 2 si 

Compute Mz and its minimal polynomial (73(2) = 121 — 132z — 36z^ + 36z'^ + 12z*. Then 
the root bound of 173(2:) is 5. We have = 10. We can set 

= i < = i 
''^ 2 - r3 2- 

Let t = X + siy + S1S2Z. Compute Mt and its minimal polynomial 

T:i{t) = 53294617 - 309903360 1 + 11884800 + 94464000 + 30720000 1'^. 



Computing its complex roots, we have 



i?3 = [-2.30803737442857- 1.39091697997219 i, -2.30803737442857 
+ 1.390916979972191, 0.174867014226204, 1.36620773463121]. 

We use i?i [i] to represent the i-tli element of Ri . R2 [i] , R3 [i] are similarly defined. Sinee 
i?2[l] - i?i[l] -0.021135190+ 0.1217738401 and the absolute values of its real part and 
imaginary part are lese than 1/2, iMlh_iMll) jg ^ j-Qot of 'P fl Q[a;,y]. But i?2[2] — 

= -0.021135190 + 2.9486217521 and its imaginary part is larger than 1/2. Then i?2[2] 
does not correspond to i?i[l]. i?3[l] -i?2[l] -0.066087950 + 0.0225069761 and the absolute 
values of its real part and imaginary part are lese than 1/4, so 

R2[l]-Ri[l] R3[l]-R2[l] ^ 

Sl ' S1S2 

= (-2.22081423399575- 1.535197796461521, -0.42270380+ 2.435476801, 
-2.64351800 + 0.900279041) 



is a root of P = 0. In a similar way, we can find all other complex roots of P ~ 0. And 
from Theorem 16, we can set ei = -^e, £2 = £3 — where e is the given precision. So if 
we refine the roots of Ti{t) = to five digits, we can obtain the roots oi V = with three 
digits. 

We also obtain an LUR for V as follows: 

[[ri(t),T2(t),r3(i)],[si,S2],[dl,d2]]. 

The roots of P = are: 

[(a,20(/3-a),40(7-/3))|Ti(a) -0,r2(/3) =0,r3(7) = 0,|/3-a| < 1/2,|7-/?| < 1/4]. 

Assuming that the final precision for the real roots of the system is e = 1/2^" and 
isolating the real roots of Ti{t) = with precision ei = ^e, £2 = £3 = ^£, respectively, we 
can obtain the following two real roots of V — with the given precision: 

,5519 345, ,4835 38695, ,20715 20725, ,44479 88959, ,-10985 -5485, ,16745 16755, 

X — — , — — X — — , — , — — , — — - X — — , — — X 



'65536 4096 4096 32768 '16384 16384' '32768 65536 ' 32768 16384 ' +6384 16384' 
In the next example, we will compare our method with RUR ([20h . 



Example 20. Consider the following example from paper (|20f). V := [ZAuz — — — 
u'^z^ — 13, 24: yz — y'^ — z^ — y'^z^ — 13, 24 uy ~ v? — y'^ — u'^y^ — 13]. The coordinate order 
is {u,y,z). 

The RUR is as follows. 

f, . ^ g{u,x) g{y,x) 9{z,x) 

j{x) — U, u ~ — -, y — — -, z = — -, 

g[l,x} g(i,x) 

where 



fix) = x^^ - 5656 + 12508972 x^^ - 14213402440 x^° + 9020869309270 

-3216081009505000 x^ + 606833014754230732 x'' 

-51316296630855044152 + 1068130551224672624689, 
g(l, x) = x^^ - 4949 x^^ + 9381729 - 8883376525 x^ + 4510434654635 

-1206030378564375 x^ + 151708253688557683 x^ - 6414537078856880519 x, 
g{u, x) = 116 x" - 483592 x^^ + 784226868 x^° - 634062241592 x® 

+270086313707548 x^ - 58355579408017944 x" + 5520988105236180668 x^ 

-131448117382500870952, 
g{y, x) = 86 x^"* - 418870 x^^ + 759804846 x^° - 670485664238 x^ + 307445009725282 x^ 

-71012402366579778 x'* + 7099657810552674458 x^ - 168190996202566563226, 
g{z, x) = 71 x^^ - 355135 x^^ + 673508751 x^° - 633214359791 x^ + 314815356659869 x^ 

-79677638700441717 x^ + 8618491509948092045 x^ - 205956089289536014429. 

An LUR of V is as follows: 

[[Ti(t),T2{t),Ta{t)], [si,S2], [di,d2]] = [[Ti(t),T2{t),Ta{t)], [1/200, 1/15], [0.2237374734, 2.146554200]], 
where 



Ti (t) = 169 - 1820 + 2622 1* ~ UOf +t^, 

Tiit) = 12034552627604020308981441166197 - 133523438810776274535699687120000 

+334257305564156882138712000000000 - 256456971612085383936000000000000 f 
+23629005541670400000000000000000 f - 665288908800000000000000000000 
+4096000000000000000000000000 t ^^ 

Tsit) = 398658124842757922827990174525891734024598098970801 

-5057045016775809265742737650285696238919118781687500 t'^ 
+18306568462902747682078658662680830721818866699218750 t* 
-26971016274307991838575084944533427932357788085937500 f 
+15563591910271113423505114668403939783573150634765625 
-1936419155067693199961145026385784149169921875000000 4^° 
+94190634217706926258139312267303466796875000000000 t^^ 
- 1851048158439662307500839233398437500000000000000 t^* 
+10022595757618546485900878906250000000000000000 



The roots of Pare: {{u,y,z) = (a, 200(/3-a), 3000(7-/3))|Ti(a) = 0,T2{f3) = 0,T3(7) = 
0, 1,3 - a| < 0.2237374734, 17- ,3| < 0.01073277100}. 



5. Conclusion 



Wc give a new representation, LUR, for the roots of a zero-dimensional polynomial system 
V and propose an algorithm to isolate the roots of V under a given precision e. For the 
LUR, the roots of the system are represented as a linear combination of the roots of some 
univariate polynomial equations. The main advantage of LUR is that precision control of 
the roots of the system is much clearer. 

The main drawback of the LUR method is that when the parameters Si becomes very 
small, the coefficients of {t) could become very large, which will slow down the algorithm. 
To improve the efficiency of the LUR algorithm is our future work. 
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